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2 Abstract. 

3 Statistical analysis of occurrence rate of magnetic storms induced by dif- 

4 ferent types of interplanetary drivers is made on the basis of OMNI data for 

5 period 1976-2000. Using our catalog of large scale types of solar wind streams 
e we study storms induced by interplanetary coronal mass ejections (ICME) 

7 (separately magnetic clouds (MC) and Ejecta) and both types of compressed 

a regions: corotating interaction regions (CIR) and Sheaths. For these types 

9 of drivers we calculate integral probabilities of storms with minimum Dst < 

10 -50, -70, -100, -150, and -200 nT. The highest probability in this interval of 

11 Dst is observed for MC, probabilities for other drivers are 3-10 times lower 

12 than for MC. Extrapolation of obtained results to extreme storms shows that 

13 such a magnetic storm as Carrington storm in 1859 with Dst = -1760 nT 
« is observed on the Earth with frequency 1 event during ~ 500 year. 
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1. Introduction 

15 One of the main problems of solar-terrestrial physics and Space Weather is study and 

16 prediction of magnetic storms including extreme ones (see, e.g., recent papers and reviews 
i? by Echer et al. [2011]; Podladchikova and Petrukovich [2012]; Yakovchouk et al. [2012]; 
is Yermolaev et al. [2012] and references therein). Useful technique for problem solution is 
is calculation of frequency (occurrence rate) distributions of events based on observations 

20 in form dN = F(x)dx, where dN is the number of events recorded with the parameter 

21 x of interest between x and x + dx, and F(x) is a frequency distributions (see, e.g., 

22 recent review and paper by Crosby [2011]; Gorobets and Messerotti [2012] and references 

23 therein). The Dst index (or its proxies) calculated from the observations of the horizontal 

24 magnetic field at four low- to mid-latitude ground stations is used to indicate value of 

25 magnetic storms [Sugiura, 1964; Sugiura and Kamei, 1991]. Recently the geomagnetic 

26 storm peak intensity distribution has been obtained for period 1957 to 2008 and it has 

27 been described by an exponential form, F(Dst p ) = 1.2exp(—Dst p /34:) where F is the 

28 probability of geomagnetic storm occurrence with a given value Dst p [Echer et al, 2011]. 

29 This formula gives too small probability for extreme storms (for instance, for storms such 
so as the 1859 Carrington storm with Dst = —1760 nT [Tsurutani et a/., 2003] probability 

31 is equal to 1 event per 10 23 magnetic storms [Echer et a/., 2011] ). 

32 To predict time and value of magnetic storm it is necessary to additionally find relation- 

33 ships between interplanetary and magnetospheric parameters/events. Alves et al. [2011] 

34 analyzed interplanetary structures (magnetic clouds with sheath before them - MC, coro- 

35 tating interaction regions - CIR, and shocked plasma) and distributions (histograms and 
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se fitting lines) for geomagnetic indices for each interplanetary structure (the geomagnetic 

37 indices Kp(ap), AE, and Dst are peak values within 2 days after the interplanetary struc- 

38 ture has passed near- Earth orbit). The main disadvantages and limitations of this work 
so are following. 

« 1. MC, Ejecta and Sheath events are not selected and analyzed together as "MC events" 

41 though it is well known that these types of interplanetary drivers result in different reaction 

42 of magnetosphere [Tsurutani and Gonzalez, 1997; Gonzalez et al, 1999; Yermolaev et al, 

43 2005; Wu and Lepping, 2007; Pulkkinen et al, 2007; Yermolaev et al, 2010c; Kilpua et 

44 al, 2012] 

45 2. The used technique of relating of interplanetary driver to magnetic storm is ambigu- 

46 ous because the delay between the reason and its consequence usually does not exceed 2 

47 hours Gonzalez and Echer [2005], and usually duration of interplanetary drivers (average 

4 8 durations 24 ± 11 h for MC, 29 ± 5 h for Ejecta, 16 ± 3 h for Sheath before Ejecta, 9 ± 
« 5 h for Sheath before MC, and 20 ± 4 h for CIR [Yermolaev et al, 2007, 2010a]) is less 

50 than 2 days. 

51 3. Used approximations at large indices (see Table 4 in paper by Alves et al. [2011]) 

52 give constant values 0.06, 0.22 and 0.045 for MC, shock and CIR, respectively, and cannot 

53 be used for analysis of extreme storms. 

54 In our previous paper [Yermolaev et al, 2012] we analyzed 798 geomagnetic storms with 

55 Dst < —50 nT and five various types of solar wind streams as their interplanetary sources: 
se corotating interaction regions (CIR), interplanetary coronal mass ejection (ICME) includ- 

57 ing magnetic clouds (MC) and ejecta, and a compression region sheath before both types 

58 of ICME (SHmc and SHej, respectively) and calculated a probability with which a se- 
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59 lected phenomenon can cause a magnetic storm, i.e., the ratio between the number of 

60 events Kj of a chosen stream type j (MC, CIR etc.) resulting in a magnetic storm with 

61 Dst < Dst and the total number of this type events Nf Pj = Kj/Nj. In contrast with 

62 papers by Alves et al. [2011] and Echer et al. [2011] we calculated an integral probability 

63 summing probabilities from — oo up to Dst , i.e. PjiDsto) = f?£° Fj(Dst)d(Dst). In 

64 previous paper we calculated P(— 50) for MC, Ejecta, Sheath and CIR. In this paper we 

6 5 calculate integral probabilities for \Dst\ = 50, 70, 100, 150, 200 nT and extrapolate them 

66 to stronger, extreme storms with \Dst\ = 500, 1000, 1700 nT. 

2. Methods 

We use our catalog of large-scale solar wind phenomena in 1976-2000 (see 

68 ftp://www.iki.rssi.ru/pub/omni [Yermolaev et al, 2009] obtained on the basis of OMNI 

69 database (see http://omniweb.gsfc.nasa.gov [King and Papitashvili, 2004] and data on 

70 Dst index (see http://wdc.kugi.kyoto-u.ac.jp/index.html). The technique of determina- 
?! tion of connection between magnetic storms and their interplanetary drivers consists in the 

72 following. If the minimum of Dst index lies in an interval of a type of solar wind streams 

73 or is observed within 1-2 hours after it we believe that the given storm has been generated 

74 by the given type of streams [Yermolaev et al, 2010]. During 1976-200 there were 798 

75 magnetic storms with minimum Dst < -50 nT but the source of 334 magnetic storms 

76 (i.e., 42 % of 798 storms) occurred to be undetermined, mainly because of the absence of 

77 a complete set of measurements for separate time intervals in the OMNI database. 

78 We use indicated data to calculate following parameters: (1) integral probability 

79 Pj(Dst ) and (2) waiting time Tj(Dst ). Pj(Dst ) = Kj(Dst )/Nj, where Nj is total num- 

80 ber of solar wind events of type j observed during period 1976-2000, Kj(Dst ) is number 
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si of solar wind events of type j induced magnetic storms with minimum Dst < Dsto, i.e. 

82 Pj(Dst ) is function of threshold Dst . In this work we calculate integral probability for 

83 values Dst = -50, -70, -100, -150, and -200 nT and for 4 types of interplanetary drivers 

84 ( CIR, Sheath, MC, Ejecta as well as ME=MC+Ejecta) on the basis of experimental 
as data. These experimental points are approximated by different functions and obtained 
ae approximations allows us to extrapolate probabilities to extreme magnetic storms with 
a? high values Dst = -500, -1000, and -1700 nT for which experimental data are too scare 
as to calculate real probabilities. 

as An important parameter of magnetic storms is a waiting time which is average period 

90 between consecutive observations of magnetic storms stronger Dsto generated by solar 

91 wind type j, i.e. Tj(Dst ) = [Pj(Dst )Nj/(t * r)] -1 , where t = 26 years (duration of 

92 observations from 1976 up to 2000) and r = 0.58 (ratio of durations of data existence to 

93 total duration of observation). To estimate waiting times Tj(Dst ) for Dst = -500, -1000, 

94 and -1700 nT, we use approximating functions pf integral probabilities at these values of 

95 Dst index. 

3. Results 

96 Total numbers of different interplanetary drivers Nj, numbers of drivers result- 

97 ing in magnetic storms Kj(Dst ) and integral probabilities Pj(Dst ) for | .Dsto I = 

98 50, 70, 100, 150, 200 nT are presented in Table 1. Data in last column with i^(-200) = 

99 3 and 5 for CIR and Ejecta were excluded from consideration of Pj because of too low 

100 statistics. Obtained data on Pj are presented in Figure 1. 

101 These data shows that the highest probability is observed for magnetic clouds. At 

102 Dsto index increase from 50 to the 200 nT the probability fell almost in 10 times. The 
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103 probability for other types is less in 3-4 times at |-Dst | = 50 and 5-15 times at 200 nT, i.e. 

104 the probability for magnetic clouds decreases more slowly with increasing index than for 
us other drives. It is important to note that ratio of probabilities of MC and Ejecta changes 
we from 5 at \Dst \ = 50 nT upto 10 at \Dst \ = 150 nT and probability of ICME (indicated 

107 as ME = MC + Ejecta) is close to probability of Ejecta, and this fact is serious argument 

108 to analyze MC and Ejecta separately. 

4. Discussion 

109 Obtained results allow us to calculate probabilities and waiting times of magnetic storms 
no in Dst range from -50 down to - 200 nT for various interplanetary drivers and to use these 
in parameters for forecasting of magnetic storms. Unfortunately, there is no enough number 
n2 of events to calculate these parameters for larger storms. Nevertheless obtained data can 
us be approximated by functions and these approximations can be used for extrapolation of 
n4 probabilities in range of very large storms. 

us In accordance with paper by Crosby [2011] various events on the Sun and the Earth have 

He similar shape of distributions: a top at small values, then a gentle slop and power-law 

n7 tail at large values. If suggestion on power-law distribution at given values of Dst is true, 

us obtained data can be approximated by power-law functions and these approximations 

n9 were made using last highest 3 measured points of \Dst\ (see Table 2 and thing lines 

120 in figure 2). As seen in Fig.2, dependence of log(P) on log(\Dst\) is not linear (as it 

121 should be for power law). So we made square approximation in log- log scales (see Table 

122 2 and thick lines in figure 2). The square-law dependencies are too steeply decreasing 

123 and do not agree conclusion by Crosby [2011] about power-low tail at large values. We 

124 suggest that dependencies do not approach power-law part of distribution in range of 



DRAFT 



November 20, 2012, 2:17am 



DRAFT 



X - 8 YERMOLAEV ET AL.: EXTREME MAGNETIC STORMS 

125 |-Dsto| = 50,70,100,150,200. So, we take power law with fixed index -2.5 for MC (see 

126 Table 2 and dashed line in figure 2) and believe that this line is the best approximation. 

127 The probabilities for another drivers in these range of Dst have been approximated by 

128 power law with fixed index -2.5 (not shown) and obtained lines look better than power-law 

129 fittings on last 3 points and square-law fittings. 

uo Probabilities obtained by extrapolations of fittings are used for calculation of waiting 

131 times for extreme storms with minimum Dst < -500, -1000, and -1700 nT (see Table 2). 

132 The data spread is very large but it is possible to make some rough estimations. The least 

133 waiting times (e.g. the highest probability) is observed for MC. Other drivers have larger 

134 waiting times (lower probabilities) and their contribution can be neglected. Thus, the 

135 most probable estimations of waiting times for extreme magnetic storms with minimum 

136 Dst < -500, -1000, and -1700 nT are ~20, ~100 and ~500 years (with accuracy of factors 

137 ~1.5, ~2, and ~3), respectively. It should be noted that we use 26 year observations to 

138 estimate waiting time of 100 and 500 years, therefore accuracy of estimates is very low. 

139 Nevertheless the obtained estimates appear more realistic than ones were obtained by 
«o Echer et al. [2011]. In order essentially to increase accuracy and reliability of estimates 
mi by means of this technique it is necessary to increase duration of an interval of observation 

142 several times. 

5. Conclusions 

143 On the basis of OMNI dataset during 1976-2000 we classified different types of solar 

144 wind events (magnetic cloud MC, Ejecta, CIR, and Sheath), found magnetic storms cor- 

145 responding to these types of interplanetary drivers and calculated integral probabilities of 
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«6 generation of magnetic storms with minimum Dst < -50, -70, -100, -150 and -200 nT by 

147 these types of solar wind drivers. Obtained data allow one to make following conclusions. 

148 1. Probabilities of storm generation by all types of drivers decrease with decreasing 

149 minimum Dst index, and rate of probability decrease increases with index decreasing. 

150 2. Probability for MC is the highest one, and rate of its probability decrease is the 

151 smallest. 

152 3. Probabilities of CIR-, Ejecta- and Sheath-induced storms with Dst < -50 nT are 

153 less in 3-4 times than for MC, for storms with Dst < -200 nT their probabilities are less 

154 in 5-15 times. 

155 Interpolation of data in Dst region from -50 down to -200 nT and than extrapolation 

156 for stronger, extreme storms with Dst < -500, -1000, and -1700 nT allow us to suggest 

157 the following. 

158 1. The most probable estimations of waiting times for extreme magnetic storms with 

159 minimum Dst < -500, -1000, and -1700 nT are ~20, ~100 and ~500 years (with accuracy 
wo of factors ~1.5, ~2, and ~3), respectively. 

mi 2. The waiting times of CIR , Ejecta- and Sheath-induced storms with Dst < -500 nT 

162 are larger in 5-10 times than for MC, their waiting times for storms with Dst < -1700 nT 

163 are larger in 10-100 times than for MC. Therefore their waiting times for extreme storms 

164 is very small relative to MC one and their contribution in extreme storm generation can 

165 be neglected. 

me Obtained experimental results and estimations are important for space weather and, in 

167 particular, for analysis and forecasting of extreme magnetic storms. 
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Figure 1. Dependence of integral probability Pj(\Dst\ > Dst ) on the value of storm for 
different types of solar wind. 
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Figure 2. Approximations of dependence of integral probability P(\Dst\ > DstO) on the value 
of storm for different types of solar wind. Vertical dashed lines show 500, 1000 and 1700 nT. 



Table 1. Probabilities Pj generation of magnetic storms with Dst < -50, -70, -100, -150 and 



-200 nT for different types of interplanetary drivers. 



Types of 
drivers 


Number, 

N S 


Dst < -50 


Dst < -70 


Dst < -100 


Dst < -150 


Dst < -200 


Kj Pj 








Kj Pi 


CIR 

Sheath 

MC 

Ejecta 

ME (MC + Ejecta) 


718 

642 

101 

1127 

1228 


120 0.167 
93 0.145 
57 0.564 
139 0.123 
196 0.159 


66 9.2 10"2 
72 0.112 
44 0.436 
89 7.9 10- 2 
133 0.108 


24 3.3 10"2 
37 5.8 10-2 
28 0.227 
45 4.0 10-2 
73 5.9 10"2 


6 8.4 10" 3 
16 2.6 10-2 
14 0.139 
12 1.1 10-2 
26 2.1 10"2 


5 

7 1.1 10-2 
7 6.9 10-2 

3 

10 8.1 10" 3 
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Table 2. Approximations of integral probabilities Pj and waiting times (average occurrence 



periods) Tj for magnetic storms with \Dst\ > 500, 1000, 1700 nT. 



Types of 


Extrapolation, 






Dst 


> 500 




Dst 


> 1000 




Dst 


> 1700 


drivers 


Y = 


log(P),X = log(\Dst\) 




P 






T 




P 




T 




P 




T 


MC 


Y = 


-1.98X + 7.87 




11.18 10 


-'2 


12.7 


2.99 10 


-3 


50 


1.05 10 


-3 


142 




Y = 


-2.5X + 4.6 




7.12 10" 


-3 


21.0 


1.26 10 


-3 


118 


3.34 10 


-4 


447 




Y = 


-4.19 + 5.14X- 


1.66X 2 


3.87 10" 


-3 


38.6 


1.93 10 


-4 


774 


1.2210 


-5 


12200 


Eject a 


Y = 


-3.26X + 5.12 




2.10 10" 


-4 


63.7 


2.20 10 


-5 


608 


3.8910 


-6 


3400 




Y = 


-8.73 + 9.91X- 


3.13X 2 


1.74 10" 


-5 


770 


7.23 10 


-8 


1.9 10 5 


4.4910- 


-10 


3.0 10 7 


ME 


Y = 


-2.85X + 4.48 




6.31 10" 


-4 


19.5 


8.77 10 


-5 


140 


1.9410 


-5 


630 




Y = 


-6.39 + 7.30X - 


2.36X 2 


1.34 10- 


-4 


91.6 


1.90 10 


-6 


6500 


3.7510 


-8 


3.3 10 5 


CIR 


Y = 


-3.42X + 5.36 




1.36 10" 


-4 


154 


1.27 10 


-5 


1650 


2.0810 


-6 


10100 




Y = 


-5.51 +7.11X- 


2.54X 2 


1.42 10- 


-5 


1480 


8.53 10 


-8 


2.5 10 5 


8.28 10 


-10 


2.5 10 7 


Sheath 


Y = 


-2.38X + 3.54 




1.30 10" 


-3 


18.1 


2.49 10 


-4 


94.3 


7.05 10 


-5 


333 




Y = 


-5.35 + 6.01X- 


1.97X 2 


3.19 10- 


-4 


73.6 


8.50 10 


-6 


2760 


3,03 10 


-7 


7.8 10 4 
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